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1.  Introduction 

The  theory  of  isotropic  random  vector  fields  was  originated  by  H.  P.  Robertson  [1] 
in  his  theory  on  isotropic  turbulence.  He  defined  the  covariance  bilinear  form  of  random 
vector  fields  which  corresponds  to  Khinchin’s  covariance  function  in  the  theory  of  sta¬ 
tionary  stochastic  processes.  Although  in  the  latter  theory  the  essential  point  was  made 
clear  in  connection  with  the  theory  of  Hilbert  space  and  that  of  Fourier  analysis,  we  have 
no  corresponding  theory  on  isotropic  random  vector  fields. 

Robertson  obtained  a  condition  necessary  for  a  bilinear  form  to  be  the  covariance  bi¬ 
linear  form  of  an  isotropic  random  vector  field.  Unfortunately  his  condition  is  not  suffi¬ 
cient;  in  fact,  he  took  into  account  only  the  invariant  property  of  the  covariance  bilinear 
form  but  not  its  positive  definite  property.  A  necessary  and  sufficient  condition  was  ob¬ 
tained  by  S.  ltd  [2].  Although  his  statement  is  complicated,  he  grasped  the  crucial  point. 
His  result  corresponds  to  Khinchin’s  spectral  representation  of  the  covariance  function 
of  stationary  stochastic  processes. 

The  purpose  of  this  paper  is  to  establish  a  general  theory  on  homogeneous  or  isotropic 
random  vector  fields,  or  more  generally  the  homogeneous  or  isotropic  random  currents 
of  de  Rham  [3].  In  section  2  we  shall  give  a  summary  of  some  known  facts  on  vector 
analysis  for  later  use.  In  section  3  we  shall  define  random  currents  and  random  meas¬ 
ures.  The  reason  we  treat  random  currents  rather  than  random  ^-vector  fields  or  £-form 
fields  is  that  we  have  no  restrictions  in  applying  differential  operators  d  and  8  to  random 
currents.  These  operators  will  elucidate  the  essential  point.  In  section  4  we  define  homo¬ 
geneous  random  currents  and  give  spectral  representations.  Here  we  shall  explain  the 
relation  between  homogeneous  random  currents  and  random  measures.  In  section  5  we 
shall  show  a  decomposition  of  a  homogeneous  random  current  into  its  irrotational  part, 
its  solenoidal  part  and  its  invariant  part.  In  the  next  section  we  shall  give  a  spectral  repre¬ 
sentation  of  the  covariance  functional  of  an  isotropic  random  current.  The  result  here 
contains  S.  ltd’s  formula  as  a  special  case.  The  spectral  measure  in  this  representation  is 
decomposed  into  three  parts  which  correspond  to  the  above  three  parts  in  the  decomposi¬ 
tion  of  a  homogeneous  random  current.  This  relation  was  not  known  to  S.  ltd.  In  [4]  we 
have  shown  that  the  Schwartz  derivative  of  the  Wiener  process  is  a  stationary  random 
distribution  which  is  not  itself  a  process.  A  similar  fact  will  be  seen  in  section  7  with  re¬ 
spect  to  the  gradient  of  P.  Levy’s  Brownian  motion  [5]  with  a  multidimensional  param¬ 
eter. 

This  work  was  supported  by  a  research  project  at  Princeton  University  sponsored  by  the  Office  of 
Ordnance  Research,  DA-36-034-ORD-1296  RD,  and  at  the  University  of  California,  Berkeley,  under 
contract  DA-04-200-ORD-355. 
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2.  ^-vector 

We  shall  here  summarize  some  known  facts  on  vector  analysis  which  will  be  used  in 
this  paper.  Let  X  =  Rn  be  the  Euclidean  w-space  and  {«,}  be  an  orthonormal  regular 
basis,  that  is, 

(2.1)  ( ®y)  =  |  if  y==l,2, 

and  the  vectors  e,-,  •••,«„,  in  this  order,  give  a  positive  orientation.  Then  any  point  xoiX 
is  expressed  uniquely  as 

(2.2)  x  =  Xid. 


The  tangent  space  Tx  and  its  dual  space  T*,  that  is,  the  space  of  differentials  at  any 
point  x  of  Xf  are  both  isomorphic  to  the  space  X  itself  by  the  following  correspondence, 


(2.3) 


dXi 


d  Xi++  ei , 


i=  1,  2,  •  •  •  ,  n  . 


Therefore,  we  may  identify  both  the  space  of  ^-vectors  at  x  and  that  of  />-forms  with  the 
space  X[pl  of  ^-vectors  with  complex  coefficients  in  X.  X  may  be  considered  as  a  real 
part  of  X[1]. 

Any  ^-vector  ap  (E  XIpl  is  expressed  uniquely  as 


(2.4)  ap—  ...  & ip, 

I<1  P  P 

where  the  coefficients  a,-,  •  •  •  i,  are  complex  numbers  and  [i]  means  that  the  summa¬ 
tion  sign  ^  refers  not  to  all  systems  of  suffixes,  but  to  those  which  satisfy  i\  <  i t  <  •  •  • 
K  ip- 

The  following  notation  will  often  be  used  in  this  paper: 
a5) 

is  equal  to  1  or  —  1  according  to  whether  {j,}  is  an  even  or  odd  permutation  of  {*„}  ,  and 
is  equal  to  0  in  all  other  cases.  We  shall  state  the  definitions  of  exterior  product  ap  Pit 
adjoint  multivector  a%  inner  product  (op,  /Sp)  and  generalized  inner  product  ap  ss  fiq, 
which  are  independent  of  the  choice  of  the  orthonormal  regular  basis, 


(2.6) 

II 

< 

s' 

*  •  ipfiif  " 

M M  [*] 

■AV- 

b  J 1 

(2.7) 

* 

aP  — 

2  “*1  ••• 

1*1  l>] 

•  •  ^ 

C12 

:  :  '  Y 

ip  j  1  *  *  *  in—p / 

(2.8)  («,,&,)  =  2  «V  <p=  *. 

1*1 

(2.9)  o,vf,  =  ( -  1)  <»-»>  (-«)  (o,^  », 

The  product  >✓  is  dual  to  ^  in  the  following  sense 

(2.10)  (a,ps-/  fiq,  yq—p)  =  (fiq)  O’p^'Yq—p)  • 


q^P  . 
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Now  we  shall  consider  a  />-vector  field  or  equivalently  a  p-form  field  a  (a;).  Although 
a  (a:),  in  its  proper  sense,  maps  each  point  x  of  X  respectively  to  a  ^-vector  or  />-form 
at  x,  it  is  also  considered  as  a  mapping  from  X  to  Xlp]  by  the  correspondence  (2.3). 
Therefore,  a(x)  will  be  expressible  as 


(2.11) 


«(*)  =  ^  *»,  •••*„(*) 

[*] 


where  the  coefficients  a»,  •  •  •  iP(x)  are  complex- valued  functions  of  x.  The  following  oper¬ 
ations  are  common  in  the  theory  of  ^-forms. 

Differential  operators  d  and  8 , 


(2.12) 


daP(x)  ek^eix^-  -^eip, 


*.  !*] 


(2.13)  8ap  (x)  =  (-D^+Udap  *. 

Inner  product , 

(2.14)  (  ap,  fipy  =  f  n  (aP>  A>)  dx!  •  •  •  dxn. 


3.  Random  current  and  random  measure 

In  this  paper  we  shall  treat  only  complex-valued  random  variables  with  mean  0  and 
finite  variance.  The  totality  of  such  random  variables  constitutes  a  Hilbert  space  which 
we  shall  denote  by  H.  The  inner  product  of  two  elements  in  H  is  the  covariance  between 
them.  We  shall  always  refer  to  the  strong  topology  in  H. 

A  random  current  is  an  H-valued  current  of  de  Rham.  Let  2>P  be  the  set  of  all 
Cm  ^-vector  fields  with  a  compact  carrier.  is  nothing  but  the  Schwartz  space. 
SbP  is  isomorphic  to  the  nCpth  power  of  ^0,  so  that  it  is  a  linear  topological  space.  A 
random  current  of  the  pth  degree  is  a  continuous  linear  mapping  from  into  H.  A 
continuous  random  ^-vector  field  Up(x)  induces  a  random  current  of  the  pth  degree  as 
follows, 

(3.1)  Up(<t>n-p )  =  fn(Up^4>n-p)  *  dx  1  •••  dxn. 

J  R 

A  sequence  of  random  currents  { U(pm) }  is  said  to  converge  to  Up  if 

(3.2)  U “  «.»-,)  Up  (*,-„) ,  €  3>n-p  . 

The  operations  on  random  currents  are  defined  in  the  same  way  as  in  the  case  of  the  cur¬ 
rents  of  de  Rham, 


(3.3) 

u;  i+j  = 

(“ 

1) 

p(n—p) 

Up  Wp) , 

(3.4) 

11 

k 

1 

(- 

1) 

p+1  u 

p  (d^>n— p— l)  > 

(3.5) 

8U— 

(“ 

1) 

np+n+1 

( dU'p) *, 

(3.6) 

aq  ( x ) 

^  Up  (fin—p—q)  = 

(“ 

1) 

M  Up 

(  0,q  /N  (f)n—p — g)  , 

(3.7) 

*q  (*) 

Up  ($„— p+a) 

(“ 

1) 

q(p—q) 

Up  <f>n— p+g)  • 

For  a  random  current  Up  we  shall  define  (Up,  ap)  as  a  random  Schwartz  distribution, 
(3.8)  ( Up,op)(<t> )  =  C/p(^-4). 
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A  random  Schwartz  distribution  M(<p)  is  called  a  random  measure  with  respect  to  a 
measure  m  if  we  have,  for  any  pair  <p,  \p  6 

(3.9)  [M  (*) ,  M  (*)  ]  H  =  fj>  (x)  TOO  «  (<**) , 

J  R 


where  [M(4>),  M(\p)]jj  denotes  the  inner  product  in  H.  Putting  M(E)  —  M(xs),  where 
Xb  is  the  characteristic  function  of  the  set  E,  we  get  a  random  set  function  which  is 
additive  in  E.  Further,  we  have 

(3.10)  [M  (E) ,  M  (Er) )  H  —  m  {EE')  ,  M  (<f>)  =  f  n<t>  (*)  M  ( dx )  . 

J  R 


Particularly,  if  m  satisfies 

(3.11) 


m(dx)  ^  , 
(1+ !*!*)*<  +" 


that  is,  m  is  slowly  increasing,  then  M (<f>)  is  also  said  to  be  slowly  increasing.  Under  this 
condition,  any  rapidly  decreasing  function  <f>  in  Schwartz  sense  [6]  belongs  to  L2(Rn,  m), 
so  that  M(4> )  can  be  defined. 

A  random  current  Mp(4>n-p)  is  called  a  random  measure  (of  the  p\h  degree)  if  there 
exists  a  complex-valued  locally  finite  measure  m(dx;  ap,  bp)  for  every  pair  (ap,  bp)  such 
that  we  have 

(3.1 2)  [  (MP)  ap)  (<#,)  ,  (Mp,  K)  =  U  (*)  TW  m  ( dx ;  a„  bp) . 

If  m(dx;  aPap)  is  a  slowly  increasing  measure  for  every  ap,  then  Mp  is  also  said  to  be 
slowly  increasing. 


4.  Homogeneous  random  current:  spectral  representation 

A  translation  ta:  x— »  x  +  h  induces  a  translation  cfh  of  p-vector  fields  in  the  follow¬ 
ing  usual  way, 

(4.1)  (.0h  •  <t>p)  ^  d Th  [<f>p  (x h)  ]  . 

By  the  identification  (2.3),  we  can  easily  see  that  drn  is  just  the  identity  mapping.  Thus 
we  have 

(4.2)  (<ta-£p)  (x)  =  4>p(x  +  h)  . 

Let  Up  be  a  random  current.  Then  we  shall  define  <ta Up  by 

(4.3)  (<rhUp)  (4>n-p)  =  Up  (o^Vn-p)  =  Up  (<7-h<t>n-p)  . 

The  functional  p(<f>p,  pp)  =  [Up(<f>*),  Up(4%)]h  is  called  the  covariance  functional  of  Up. 
The  function  defined  by 

(4.4)  p  (<f>,  ip)  ap>  bj)  =  [(Up,  flp)  (<p)  ,  (Up,  b p)  (\p)  ]  jj 

is  called  the  covariance  bilinear  form  of  Up  and  is  a  generalization  of  the  form  Robertson 
used  in  his  theory  of  turbulence.  If  the  covariance  functional  of  crh Up  is  independent  of 
h,  then  Up  is  said  to  be  homogeneous.  If  this  condition  is  satisfied,  the  covariance  bilinear 
form  of  o- a  Up  is  independent  of  h.  The  converse  is  also  true. 

As  in  the  theory  of  stationary  stochastic  processes  we  have  the  following  theorem 
of  spectral  representation. 
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Theorem  4.1.  The  covariance  bilinear  form  of  any  homogeneous  random  current  is 
written  as 

(4.5)  p  (<f>,  yp;  ap,  bp)  =  fR„  #<P  (y)  (y)  m  ( dy ;  ap,  bp) 

where  m(  A;  &p)  is  a  positive  definite  bilinear  form  in  (bp,  ap)  and  m(  A;  ap,  ap)  is  a  slow¬ 
ly  increasing  nonnegative  measure.  Conversely  the  p  defined  by  (4.5)  is  the  covariance  bilinear 
form  of  a  homogeneous  current. 

Remark.  #4>  is  the  Fourier  transform  of  <f>,  that  is, 

(4.6)  (4fy)  (y)  =  f  ne~i2T(x>  y)<f>(y)  dy . 

For  a  ^-vector  field  4>.(y),  we  define 

(4.7)  (<?■*.„)  (y)  •  ip(y)  «fj). 

i»i 

Theorem  4.2.  .4  homogeneous  random  current  Up  is  the  Fourier  transform  of  a  slowly 
increasing  random  measure  Mp,  which  is  called  the  spectral  measure  of  Up,  that  is, 

(4.8)  Up  (<t>n-p)  =MP  =  MP  {4<bn-P) . 

We  can  prove  this  theorem  easily  by  remarking  that  ( Up ,  aP)(<p)  is  a  stationary  ran¬ 
dom  distribution  [4]  with  a  multidimensional  parameter. 


5.  Homogeneous  random  current:  canonical  decomposition 

In  this  section  we  shall  discuss  a  decomposition  of  a  homogeneous  random  current  Up 
into  an  irrotational  current,  a  solenoidal  current  and  an  invariant  current.  This  de¬ 
composition  will  be  called  the  canonical  decomposition  of  Up.  We  shall  start  from  the 
spectral  representation.  By  using  this,  we  can  show  the  existence  of  the  limit 

(5.1)  A/Up(<t>n-p)  =  lim  f  •••  f  <rhUP(<l>n-p)  dhi dhn. 

A—* co  AnJ-A  J-a 


We  say  that  ZJP  is  invariant  or  unbiased  according  to  whether  MUP  =  Up  or  0.  An 
unbiased  homogeneous  random  current  is  called  irrotational  or  solenoidal  according  to 
whether  dTJp  —  0  or  8UP  =  0.  We  set 

(5.2)  M°p(. A)  —Mp  (An  {0}),  M$(  A)  =  MP  (A  —  { 0 } ) . 


Then  we  have  MU p  =  %MP. 

Theorem  5.1. 

(a)  A  homogeneous  random  current  U  is  irrotational  if  and  only  if  Mp  =  0  and 
y  Mp(dy)  =  0. 

(b)  U  is  solenoidal  if  and  only  if  Mp  =  0  and  y  \/  Mp(dy)  —  0. 

The  essential  point  of  the  proof  is  as  follows. 

(5.3)  d  Up  (<t>n-p- 1)  =  (  -  1)  *+1  Up  Wn-p-x)  =  (  -  1)  >+1Mp  (Qdtn-p-!) 

=  (  -  1) p+l  Mp  (i2vy  4^»-p-i)  . 


Now  we  shall  introduce  two  random  measures  Mp  and  Mp  by 


(5.4) 


M‘,(dy)  = 


y  / s 


[y^s  M%(dy)  1 

\y  I* 


[y^M^(dy)  1 

\y\ 


s 
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Theorem  5.2.  Up  =  Up  =  #MP  and  U%  =  4}MP  =  MUP  are  respectively 

irrotational,  solenoidal  and  invariant ,  and  we  have 

(5.5)  Up=  U'p+U‘p+U°p. 


By  using  the  following  identity 

t'J.o;  ap  ap  i  ap—  |y|2  '  \  y  \  2  ’  ®p>  ®p  ^  » 

we  can  prove  the  theorem. 

Now  we  shall  define  Wp(<f>n-P)  by  the  following  procedure, 


(5.7) 


G(x,y)  = 


”  g-2w-(*.  V)  — 

1  +  2ir  (x,  y) 

if 

|  y  |  <  1 , 

e-2ir(x,  v) 

y\ 2 

1  I?  I2 

if 

\y\  ^ i ; 

(5.8)  G  (4>p,  y)  =2  LfG(s,  y)  <*>,,•••  ip(x)  dx]  e^-  •  •  /\  eip; 

[i] 

(5.9)  W„  (<(..-„)  =J^KIG  (*.-2,  y)  ] . 

Theorem  5.3. 

(5.10)  d8Wp=U'p,  8dWp  —  Up . 

Therefore  the  canonical  decomposition  is  written  as 

(5.11)  Up=  d8WP+  8dWp+  U°p  . 


6.  Isotropic  random  current 

To  begin  with,  we  shall  introduce  some  preliminary  notation.  Let  G  be  the  whole 
group  of  orthogonal  transformations  (with  determinant  ±1)  in  X  —  Rn.  We  can  carry 
out  the  same  procedure  for  g  6  G  as  we  did  for  the  translation  rh  in  section  4.  Corre¬ 
sponding  to  (4.1)  we  have 

(6.1)  (<To<t>P)(x)  =  dr0  [<t>p(g'x)  1  . 

Although  dr0  here  is  not  the  identity  mapping,  we  can  easily  see  that  dr0  —  f1.  There¬ 
fore,  we  have 

(6.2)  (<W„)  (*)  =  g"1  •$,(*•*)  , 

where  is  defined  by 

(6.3)  g~l<t>p= 

[*] 

Generalizing  this  transformation  <r0 ,  we  define  <ra  on  a  current  up, 

(6.4)  ( <T oWp)  (0n— p)  =  Wp  ( O'ff-10n— p)  =  Wp  [ g $n— p  (§  1#)  1  • 

cr0  is  clearly  commutative  with  the  differential  operators  d  and  8. 

Now  we  shall  define  an  isotropic  random  current.  A  homogeneous  random  current  is 
said  to  be  isotropic,  if  the  covariance  functional  of  <r0Up  is  independent  of  g  £  G.  Even 


ISOTROPIC  RANDOM  CURRENT  131 

if  we  replace  the  covariance  functional  by  the  covariance  bilinear  form  in  the  above 
statement,  we  shall  obtain  an  equivalent  definition. 

Let  Up  be  an  isotropic  random  current.  Since  Up  is  homogeneous,  the  covariance  bi¬ 
linear  form  is  written  as 

(6.5)  p  (<t>,  p;  ap,  bp)  =  f#<f>4pm  ( dx ;  ap,  bp) 
by  theorem  4.1.  Since  Up  is  isotropic,  we  obtain 

(6.6)  tn(g-dx;  g-ap,  g-bp)  =  m(dx\  ap,  bp) . 

By  making  use  of  this  property  we  obtain  the  following  theorem. 

Theorem  6.1.  In  an  isotropic  turbulence  the  m(dx\  ap,  bp)  in  (6.6)  is  expressible  as 

(6.7)  m  (dx;  ap,  bp)  =  (0v  b^d^ap)  dd  F,  (dr) 

+  (0^  bp,d/\ap)  ddF,(dr)  +  (bp,  ap)  F0  (dx) , 

where  r  —  |  x  | ,  6  =  x/r ,  x  7*  0  ,dd  is  the  surface  element  of  the  unit  sphere ,  Fi  and  F2  are 
both  slowly  increasing  nonnegative  measures  on  (0,  °°),  and  F 0  is  a  nonnegative  measure 
such  that  F0(A)  =  0  for  0  $  A. 

Conversely ,  p(<t>,  p)  ap,  bp),  if  determined  by  m(dx\  ap,  bp)  of  the  form  (6.6),  is  the  co- 
variance  bilinear  form  of  a  certain  isotropic  random  current. 

Remark.  According  to  whether  p  =  0  or  p  =  n,  the  first  or  the  second  term  in  (6.7) 
will  disappear. 

We  shall  sketch  the  proof.  First  we  consider  the  case  in  which  m(dx ;  ap,  bp)  has  a 
continuous  density /(*;  ap,  bp).  Here  f(x\  ap,  bp)  is  a  positive  definite  bilinear  form  in 
(bp,  ap)  and  is  invariant  in  the  sense  f(gx;  gap,  gbp)  =  f(x;  ap,  bp).  If  we  introduce 
F(x;  &,•••,  £p,  771,  *  *,  Vp)  =  /(*;  Si  ^  ^  h,  rn  ^  ^  Vp),  for  real  1-vectors  rjjy 

we  obtain  a  function  which  is  linear  in  £»  and  in  rjj  and  skew  symmetric  in  each  of  {£,•} 
and  {j;,}.  Since  F  is  invariant  under  orthogonal  transformations,  it  is  a  function  of 
(£»,  &)»  (ti,  Vi),  (Vi,  Vi),  (*,  ti),  (x,  Vi)  and  (x,  x).  Using  these  properties  we  can  prove 
that /  can  be  written  in  a  form  similar  to  (6.7). 

We  can  discuss  the  case  of  general  measure  by  approximating  it  by  measures  with 
continuous  density  in  Helly’s  sense. 

The  following  theorem  shows  that  the  decomposition  in  (6.6)  corresponds  to  the 
canonical  decomposition  (5.5). 

Theorem  6.2.  In  the  case  of  an  isotropic  random  current  Up,  Up,  Up  and  UP  are  all  iso¬ 
tropic  and  orthogonal  to  each  other  in  H.  The  m-measures  corresponding  to  these  three 
parts  are  respectively  the  three  parts  in  the  decomposition  (6.7). 

7.  The  case  of  p  =  1 

In  the  case  of  p  =  1,  the  decomposition  (6.7)  becomes  simpler, 

(7.1)  m(dx;  a, /3)  =  aMj3y  [  6p6vd  6  Fj  (d  r)  +  ( —  6pQy)  Ft  (d r)  H~ F0  (dx)  j  , 

M.  v 

where  a  =  apeM  and  /3  = 

In  the  case  of  isotropic  1-vector  fields  the  measures  F,  and  F,  will  be  bounded. 

(a)  Robertson's  isotropic  turbulence.  This  is  an  isotropic  random  1-vector  field  in  3- 
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space.  Our  decomposition  (7.1)  is  essentially  the  same  as  S.  ltd’s  formula  (2)  but  is  some¬ 
what  simpler. 

(b)  The  gradient  of  P.  Levy's  Brownian  motion  with  a  multidimensional  parameter. 
Let  B(x)  be  Levy’s  Brownian  motion  with  an  n-dimensional  parameter  [5].  Then  we  have 

(7.2)  [B{x),B{y)]H=\{\x\  +  \y\-\x-y\). 

Therefore,  by  simple  computation  we  can  obtain  the  following  expression  of  the  covari¬ 
ance  bilinear  form  of  the  gradient  dB  of  B, 

(7.3)  p  (<j>,  p,  a,  j8)  =  fn#<t>  (x)  (x)  m  ( dx ;  a,  0)  , 

J  R 

where 

(7.4)  m(dx;  a,  /S)  =C„^  aMj 8,dM6,ddr~(n~l)dr , 

M.  » 

where  r  —  \x\,d  —  x/r,  and  Cn  is  a  positive  constant.  This  is  the  special  case  of  (7.1) 
in  which 

(7.5)  Fi  (d  r)  =  Cn  •  r~(n-1)  d  r  ,  F.  =  F0=0. 

Since  Cn*r~<n~^dr  is  an  unbounded  measure,  dB  is  a  random  proper  current,  that  is,  a 
random  current  that  is  not  itself  a  random  1-vector  field.  Since  F,  =  Fo  =  0  we  see  that 
dB  is  irrotational. 
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